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ABSTRACT 



Based on the accretion induced magnetic a field decay model, in which a frozen field 
and an incompressible fluid are assumed, we obtain the following results. (1) An analytic 
relation between the magnetic field and spin period, if the fastness parameter of the accre- 
tion disk is neglected. The evolutionary tracks of accreting neutron stars in the P-B dia- 
gram in our model are different from the equilibrium period lines when the influence of the 
fastness parameter is taken into account. (2) The theoretical minimum spin period of an 
accreting neutron star is max ^l.lms (f^-) R^^^hs (j^) ^ > 1.1ms (jjf^) ^ R 1 ^ 1 ^ , 
independent of the accretion rate (X-ray luminosity) but dependent on the total accretion 
mass AM. However, the minimum magnetic field depends on the accretion rate. (3) The 
magnetic field strength decreases faster with time than the period. 

Subject headings: stars: neutron - stars: rotation - stars: magnetic field -X-rays: 
stars-stars: binaries: close 



2 



1. Introduction 



The magnetic field of a neutron star has long been a complex issue and one which is 
not yet to be solved(Bhattachaya & van den Heuvel 1991; Chanmugam 1992; Phinney 
& Kulkarni 1994). On the evolution of the magnetic field of a neutron star, there is not 
yet a commonly accepted model (for a general review cf. Bhattacharya and Srinivasan 
1995). The currently popular idea seems to ascribe the field decay of a neutron star in 
an X-ray binary to the period during which accretion occurs. There is evidence that 
the magnetic fields of X-ray neutron stars and recycled pulsars are correlated with the 
duration of the mass accretion phase, or the total amount of matter accreted(Taam & 
van den Heuvel 1986; van den Heuvel et al 1986). In fact, Taam and van den Heuvel 
have already discovered a possible inverse correlation between the magnetic field and the 
estimated total mass of accreted matter for the binary X-ray sources. Later, Shibazaki 
et al. (1989) presented an assumed formula relating the decay of magnetic field with 
accretion mass,which seems to reproduce the observed field-period relations of the recycled 
pulsars quite well. 

Theoretically, for explaining the accretion induced field decay, some suggestions 
and models have been proposed(e.g. Bisnovatyi-Kogan & Komberg 1974; Romani 1990; 
Ruderman 1991a,b,c; Dinget al. 1993; Zhang et al. 1994; Urpin & Geppert 1995; Geppert 
et al 1996; Urpin & Konenkov 1997; Zhang et al 1997; Cheng & Dai 1997; Zhang 1998; 
Ruderman, Zhu & Chen 1998). Recently, van den Heuvel & Bitzaraki(1995a, 1995b), 
from the statistical analysis of 24 binary radio pulsars with nearly circular orbits and low 
mass companions, discovered a clear correlation between spin period and orbital period, 
as well as between the magnetic field and orbital period. These relations strongly suggest 
that an increase in the amount of accreted mass leads to a decay of the magnetic field, 
and a 'bottom' field strength of about 10 8 G is also implied. White & Zhang(1997) 
discovered that the spin periods of LMXBs, implied by killo-hertz X-ray QPO, constitute 
a homogenouse group with spin period of about 2 milliseconds, which has little correlation 
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with X-ray luminosity. The above two recent observational statistics seem to place some 
constraints on the construction of a theoretical model of accretion induced magnetic field 
decay. 

In this paper, we study the evolution of a magnetic field and spin period of accreting 
neutron stars according to an accretion induced magnetic field decay model(Zhang et al 
1997; Cheng & Zhang 1998), based on the idea of van den Heuvel & Bitzaraki(1995a, 
1995b) and Romani(1990) for the mechanism of accretion induced field decay. In this 
model, the accretion matter starts to be channeled onto the two polar caps by the strong 
magnetic field near the Alfven radius. Part of the accreted matter flowing towards the 
equator pushes the field lines aside and thus dilutes the polar field strength. The bottom 
field should be reached when the polar cap extends over the entire stellar surface, which 
corresponds to the Alfven radius matching the star's radius, and gives a stellar magnetic 
field of about 10 8 G. Since the spin period of an accreting X-ray neutron star depends sen- 
sitively on the magnetic field together with the influence of the fastness parameter, some 
interesting results can be obtained, which are consistent with the recent observational 
data on low mass X-ray neutron stars by White & Zhang (1997) and on the millisecond 
pulsars by van den Heuvel & Bitzaraki(1995a, 1995b). 



2. Models 



Under the assumption that the magnetic field lines of the neutron star are frozen in the 
entire crust which has homogenous average mass density, we (Cheng & Zhang 1998) have 
obtained an analytical expression for the field evolution as follows, 

D = B± (1) 

{1 - Cexp(-fg)}^ 1 } 

where C = 1 — Xq and Xq = ("§^) 2 ^ 7 , Bq is the initial magnetic field strength, M cr is 
the crustal mass, AM = Mt and Bf is the magnetic field defined by the Alfven radius 
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matching the radius of neutron star, i.e., Ra(Bj) = R, which gives, 

B f = i.Sxltf(My<*(£y'*R?'< G, (2) 

M Ed IVIq 

where M Ed = 10 18 i? 6 gs _1 is the Eddington accretion rate and R 6 is the radius of the 
neutron star in units of 10 6 cm. 

Three consistent observational conclusions can be obtained, (1) the field decay is 
inversely related to the accreted mass, (2) the bottom field strength is about 10 8 Gauss, 
and (3) the bottom field strength is proportionally related to the X-ray luminosity. On the 
basis of the above solution, we will study the spin period evolution of accreting neutron 
stars. 

To acquire the magnetic field versus period relation, we use the formula for the 
variation of the rotation due to accretion given by Ghosh & Lamb (1979, hereafter GL) 

- P = 5.0 x l^^r^R^&Bl^P^fnM syr" 1 , (3) 

where B\ 2 is the surface field in units of 10 12 G, Z45 is the moment of inertia in units of 
10 45 gcm 2 , L37 is the X-ray brightness in units of 10 37 ergs _1 and n(u s ) is a dimensionless 
function that depends primarily on the fastness parameter, 

^ = im^r 2/7 Rl 5/7 ]B% 7 p- 1 L^ 7 - (4) 

For a star rotating slowly in the same sense as the disk flow (u s <C 1), GL found that 
n{u) s ) ~ 1.4. They also found that the dimensionless function n(u s ) decreases with 
increasing u s and becoming negative for oj s > u c . A simple expression for n(u> s ) that 
agrees approximately with numerical results over the whole range of u s is, 

„(.,) = , 4 x . (5) 

GL found uj c ~ 0.35 from their model, but stressed that the actual value of this criti- 
cal fastness parameter was relatively uncertain. Subsequent work(Ghosh & Lamb 1991) 



indicates that oo c is unlikely to be less than 0.2, but it could be as large as 0.9. In the 
following subsections, we study the influence of the field decay on the evolution of the 
spin period of the neutron star in an X-ray binary. 



2.1 Spin Evolution with a Constant Fastness Parameter 



If the variation of the fastness parameter can be ignored, i.e. n(ou s ) =1, we can 
solve equation (3) analytically. First we can rewrite equation (3) as j£r = —5.0 x 
lO-^^r^Rl^I^Bf^L^z^dtsji- 1 , where B fl2 = B f /10 12 G and z = B/B f . 
From equation (1), we obtain dt = ~ ^ cr z J 7 Z _ V which gives % = 10" 4 ^(^)" 3/7 ' Rf n 
B^{ 7 2 L% 7 syr -1 , where y = z 2 ^ 7 . Using equation (2) to eliminate Bf, assuming mass 
and radius as constants and the initial period condition P = 00 , we obtain the field- 
period relation (B-P) in the following analytic form, 

L5U1H (M/M )- 1 / 2 i2 6 - 5/M / 4B (M cr /O.lM )- 1 (6) 



atanh(x) — atanh(x ) 



where x = y' 1 = {^-) 2/7 ,x = (|^) 2/7 and atanh(x) = [ln(l + x) - ln(l - x)]/2. There 
are two interesting limits for equations (1) and (6). For AM <C M cr , equations (1) and 
(6) can be approximated as 

B « Bj^Z- oc r 1 (7) 
AM 



and 



1.5 ms ( M \ 1/2 „_, /14 . / M„ x 1 



For AM 3> M cr , equations (1) and (6) can be approximated as 
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and 

We should note that the period of the neutron star cannot be shorter than the equilibrium 
spin-up line which represents the minimum period (Bhattacharya & van den Heuvel 1991) 

- 3 / 7 / \ - 5 /7 



M \ 1 ( M 



P eq = 2AmsB 9 [^} i^ 6/7 , (11) 



to which such a spin-up may proceed at the Eddington accretion rate. If we substitute 
the minimum B field in equation (11), we obtain the minimum equilibrium period, 

PZn = 1-lms (^)~ V V 14 . (12) 

However, we want to point out that equation (12) is valid only in cases where the neutron 
star has accreted a sufficient amount of matter to spin-up to that period. For a given 
amount of accreted matter, the minimum period of the neutron star is 



Pmin max 



an expression which is independent of the accretion rate (X-ray luminosity) but depends 
on the total accretion mass AM. 

We plot the field-period relation in figure 1, which shows the evolutionary track 
curves of the magnetic field and spin period (curve 1 and curve 3). Solid diamonds in 
figure 1 are observed data summarized in table 1. Initially, a small amount of mass is 
transferred, and the neutron star is spun-up from the death valley (Chen & Ruderman 
1993) where it has a long period, which causes a modest field decay and produces systems 
such as PSR0655+64 and PSR1913+16. The binaries with longer-lived accretion phases, 
e.g.LMXB, will accept sufficient mass from their companions, and yield a substantial field 
decay as in the case of millisecond pulsars such as PSR1953+10 and PSR1620+21. Our 



analytic model B-P curves, which are obtained by assuming n{ui s ) = 1, can go beyond the 
spin-up line (the equilibrium period line). This results from the fact that the influence 
of the fastness parameter is neglected, therefore the spin-up torque still exists even when 
the spin of the neutron star is faster than the Keplerian angular velocity of the accretion 
disk at the inner radius. /,From the analytical P-B relation, we find that the minimum 
period can be obtained if the magnetic field arrives the bottom field strength(cf. equation 
6). The maximum mass accreted in LMXB from the companion could be ~ l.OM (van 
den Heuvel & Bitzaraki 1995a, 1995b), the mass of the neutron star after accretion could 
reach ~ 2.4M with a radius R§ ~ 1 for realistic equations of state (cf. Table I of Cheng 
& Dai 1997) and I 45 oc M. So the minimum spin period given by equation (10) is about 
Pmin ~ 1.7 ms. Further, the more interesting thing is that, unlike the bottom field, the 
minimum period is independent of the accretion rate (X-ray luminosity) (cf. equations 
(10) and (12)). Our expression seems to be supported by the recent work of White & 
Zhang(1997). They find that the luminosities in the 10 samples of QPO LMXBs vary 
by two order of magnitude from L 36 = 1 to L 38 = 1, but the spin periods of the sources 
diffuse into a narrow region from 2.76 ms to 3.8 ms. 

2.2 Spin Evolution with a Non-constant Fastness Parameter 

Numerical solutions (curve 2 and curve 4) for equations (1), (3), (4) and (5) where 
the critical fastness parameter is set at 0.9 are plotted in figure 1. However, we find that 
the influence of the fastness parameter has little effect in the low magnetic field region. 
The main effect of the fastness parameter is to force the evolution curves back/below 
the equilibrium period line. The evolutionary curves in the P-B diagram, however, are 
not sensitive to the fastness parameter near the bottom field and/or near the minimum 
period. /,From figure 1 these curves show that the field decay time scale is longer than the 
spin-up time scale because the field reaches the bottom value first, then the neutron star 
evolves towards the minimum spin period horizontally in the P-B diagram. The numerical 
solution are consistent with our analytic expressions in equations (7) and (8). 
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However, the fastness parameter effect is important when the evolutionary track is 
close to the equilibrium period line, which means that the spin angular velocity of the star 
matches the Keplerian angular velocity at the inner edge of the accretion disk, and the 
accretion torque produced by magnetic lines immersed in the accretion disk tends to pro- 
duce the negative (spin-down) torque. This effect ensures that the evolution track cannot 
go beyond the equilibrium period line. It is interesting to note that none of the model B-P 
curves go along the equilibrium period line, which may seem confusing. However, in fact, 
the equilibrium period line only represents the final position of the evolution track in the 
P-B diagram. Therefore it is not surprising that the real evolutionary track deviates from 
the equilibrium period line. Physically, the accretion induced field decay arises from the 
contraction of the corotation radius of the magnetosphere during the accretion spin-up 
phase, and the evolution track has little chance to meet the equilibrium line if the field 
decay time scale is shorter than the spin-up time scale at the early stage of the accretion 
phase. Some X-ray sources in HMXB such as Her X-l and Vela X-I should be very close 
to the equilibrium period position if the field decay really exists in the accretion phase. 



3. Conclusion 



We have presented a simple model for the evolution of the magnetic field and spin 
period of accreting neutron stars. Analytic formulae for evolution trajectories in B-P 
are derived. The theoretical minimum period of the neutron star does not depend on 
the accretion rate but instead depends on the total amount of accreted matter, and the 
stellar parameters including the moment of inertia, stellar mass and radius, which depend 
on the equations of state. Our model results seem to be supported by the observed data 
(White & Zhang 1997). However, in this paper we have ignored the fact that the stellar 
parameters, i.e. M, R and J, are all time dependent. The exact evolution curves must 
take this factor into account. On the other hand, the minimum period should not depend 
on the details of the evolution trajectories instead it only depends on the final values of 
the stellar parameters which are equations of state dependent. Observing the minimum 
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period in LMXB may provide useful constraints on the equation of state for high density 
matter. For AM < 1M , the minimum period is longer than 1.7ms for a wide range of 
realistic equations of state (Wiringa, Fiks & Fabrocini 1988). 

We thank A.Potekhin for his useful comments and P.K.MacKeown for a critical 
reading of our manuscript. This research program is partially supported by a RGC grant 
of the Hong Kong Government, a Croucher Foundation Senior Research Fellowship and 
a NSF grant of PRC. 
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Figure 1: Magnetic field vs. spin period diagram. Solid diamonds are observed data 
summarized in table 1. Initial spin period is chosen to be 100 s. The initial magnetic 
fields B = 5 x 10 12 G and B = 1 x 10 12 G are used, and the luminosity is varied from the 
Eddington limited luminosity £38=1 to L 36 =l. The upper heavy solid line is the spin- up 
line (equilibrium period line) given in equation (11) and the lower solid line is the death 
line defined as B 12 /P 2 = 0.2. Curve 1 and curve 3 are calculated from equation (6) with 
luminosity L 38 = 1 and L 36 = 1 respectively. Curve 2 and curve 4 are our numerical 
solutions for equations (l)-(5) with luminosity L 38 =l. and L 36 =l. respectively. In all 
model evolution curves, we have used M = 1M Q and R = 10 6 cm. 
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